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1 Introduction 

Let ft be a domain is (N > 2) with a C 2 compact boundary <9f2. A function u £ W,^(Q) 
is p-harmonic if 

/ \Du\ p ~ 2 (Du : D<t>)dx = (1.1) 
Jn 

for any </> £ Cg(ri). Such functions are locally C 1,a for some a £ (0, 1). In the case p = N, 
the function u is called A^-harmonic. The Af-harmonic functions play an important role as a 
natural extension of classical harmonic functions. They also appear in the theory of bounded 
distortion mappings 8 . One of the main properties of the class of A^-harmonic functions 
is its invariance by conformal transformations of the space M. . This article is devoted to 
the study of Af-harmonic functions which admit an isolated boundary singularity. More 
precisely, let a £ dft and u £ Wf^ c (O) PI C(Q \ {a}) be a A^-harmonic function vanishing on 
dft \ {a}, then u may develop a singularity at the point a. Our goal is to show the existence 
of such singular solutions, and then to classify all the positive A^-harmonic functions with a 
boundary isolated singularity. We denote by n a the outward normal unit vector to f2 at a 
The main result we prove are presented below: 

There exists a unique positive N -harmonic function u = u\ ta in Q, vanishing on dfl \ {a} 
such that 

lim \x — a\ u(x) = — (a, n a ) (1-2) 

x — > a 

uniformly on S N - 1 nO = {a £ S N - 1 : (a, n a ) < 0}. 

The functions u\ ta plays a fundamental role in the description of all the positive singular 
A-harmonic functions since wc the next result holds 

Let u be a positive N -harmonic function in Q, vanishing on dfl \ {a}. Then there exists 
k > such that 

u = fcwi, a - (1-3) 
*To appear in Communications in Partial Differential Equations 
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When u is no longer assumed to be positive we obtain some classification results provided 
its growth is limited as shows the following 

Let u be a N -harmonic function in Q, vanishing on \ {a} and verifying 

\u\ < Mui >a , 

for some M > 0. Then there exists k £ K such that 

u = kui ia . (1.4) 

In the last section we give a process to construct p-harmonic regular functions (p > 1) 
or TV-harmonic singular functions as product of one variable functions. Starting from the 
existence of p-harmonic functions in the plane under the form u{x) = u(r,o~) = r"ui(9) 
(see 0), our method, by induction on TV, allows us to produce separable solutions of the 
spherical p-harmonic spectral equation 

// o\ (P~2)/2 \ / „\ (p-2)/2 

-div a H/3V + |V CT «| 2 J V a v\ = \ N , p [f3 2 v 2 + \V a v\ 2 j v. (1.5) 

on S^ -1 , where \n,p = P (N — 1 + (f3 — l)(p — 1)). This equation equation is naturally 
associated to the existence of p-harmonic functions under the form u(x) = \x\^ v{xj \x\). As 
a consequence, we express p-harmonic functions under the form of a product of A^-explicit 
functions of one real variable. If we represent M. N as the set of {x — (xi, x^)} where 
x\ = r sin^jv_i sin^jv-2-" sin^2 sin^i, X2 = rsin0jv-isin#jv-2"-sin#2Cos0i, xn-i — 
r sin (9^! cos 9n~2 o-nd xn — t cos 0jv-i with 9\ G [0, 27r] and 9u £ [0, tt], for k — 2, N — 1, 
then, for any integer k the function 

u(x) — (r sin^jv-i sin0 2 )^ fc ^fc(^i) (1-6) 

is p-harmonic in M. N , in which expression (3k > 1 is an algebraic number depending on k 
and LOk is a it / k-antiper iodic solutions of a completely integrable homogeneous differential 
equation. Moreover N -harmonic singular functions are also obtained under the form 

u(x) =r- /3k (sm9 N ^ 1 sm9 N ^2-sm9 2 ) l3k Lu k (9 1 ). (1.7) 

Our paper is organized as follows: 1- Introduction. 2- Construction of fundamental 
singular TV-harmonic functions. 3- The classification theorem. 4- Separable solutions of the 
p-harmonic spectral problem. 

2 Construction of fundamental singular iV-harmonic func- 
tions 

We denote by Hn the group of conformal transformations in R N . This group is generated 
by homothethics, inversion and isometries. Our first result is classical, but we repeat the 
proof for the sake on completeness. 

Proposition 2.1 Let u be a N -harmonic function in a domain G C M. N and h £ 7ijq . Then 
Uh = uo h is N -harmonic in h~ 1 (G). 
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Proof. Because for any p > 1 the class of p-harmonic functions is invariant by homothethies 
and isometries, it is sufficient to prove the result if h is the inversion Tg with center the 
origin in R N and power 1. We set y = Iq(x) and v(y) = u(x). For any i = 1, TV 

j 

Then 

| J D U | 2 ( 2 ;) = | a ;r 4 |^| 2 ( 2/ ) = | 2/ | 4 | J Dt;| 2 (y). 
If is a test function, we denote similarly tp(y) = 4>{x), thus 

(Du, D(f>) = \x\~ 4 (D«, £>V) = I2/I 4 

and 

/ \Du\ N ~ 2 (Du 7 D(j))dx = [ \y\ 2N \Dv\ N ' 2 (Dv,D^)\Dl^\ dy 
Jg J il(G) 

Because = \det{dxi/ dyj)\ — \y\ 2N , the result follows. □ 



Proposition 2.2 Let TV > 2, B — i?i(0) and a G dB. Then there exists a unique positive 
N -harmonic function U l in B which vanishes on dB \ {a} and satisfies 

U l {x) = ^— !fL(i + (i)) asx^a. (2.1) 
\x — a\ 

Proof. We first observe that the coordinates functions xi are TV-harmonic and positive in 
the half-space Hi = {x e : Xi > 0} and vanishes on dHi. Therefore, the functions 
Xi( x ) = x i/ \ x \ are a l so TV- harmonic and singular at 0. Without loss of generality we can 
assume that a is the origin of coordinates, and that B is the ball with radius 1 and center 
(—1,0, ...,0). Let u> be the point with coordinates (—2,0, ...,0). By the inversion Z^, a is 
invariant and B is transformed into the half space H±. Since xi is TV-harmonic in Hi, the 
function 

T4, \ \ x \ 2 + 2x l 
x * XioI u (i) = 2 

2 |x| 

is TV-harmonic and positive in B = {x : | cc | 2 + 2xi < 0}, vanishes on dB and is singular at 
x = 0. If we set x[ = Xi + 1, x- = X, for i = 2, TV and U l {x') = Xi ° Zt( x )> then the x' 
coordinates of a are (1, 0, 0) and 



^) = ^T7^i = lV J 4( 1 + °C 1 )) as;c ' 

2 \x' - a x' - a 



Let t/ 1 be another positive TV-harmonic function in T3 which verifies (|2.1 (I and vanishes on 
dB \ {a}. Thus, for any S > 0, (1 + 8)U i , is positive, TV-harmonic, and U l — (1 + <5)[/ 1 is 
negative near a. By the maximum principle, U % < (1 + 8)U l . Letting S — > 0, and permuting 
t/ 1 and [7 l yields C/ 1 = f/\ □ 

By performing the inversion Zq , we derive the dual result 
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Proposition 2.3 Let N > 2, G = Bf(0) and a G dB. Then there exists a unique positive 
N -harmonic function U e in G which vanishes on dB \ {a} and satisfies 

U e {x) = o(ln|x|) as|x|-K», (2.2) 

and 

U e (x) = |X| - ^ (1 + 0(1)) 8S ^a. (2.3) 
|x — a| 

Proof. The assumption (|2.2 I) implies that the function U = U e o T\ , which is iV-harmonic 
in B \ {0} verifies 

U(x) = o(ln(l/|x|)) near 0. 

By is a removable singularity and thus U can be extended as a positive iV-harmonic 
function in B which satisfies (|2.1 [) . This implies the claim. □ 

We denote by p(x) the signed distance from x to dfl. Since c?Sl is C 2 , there exists /3o > 
such that if x G M. N verifies — (3q < p(x) < 0o, there exists a unique £ x € 917 such that 
I a: — £ x \ = |p(ar) | . Furthermore, if is the outward unit vector to dfl at £ x , x = f^ — p(x)u^ :c . 
In particular £3 — p(x)f£ x and ^ + ^(x)^^ have the same orthogonal projection £ s onto <9il. 



Let I> (n) = {a; G l w : -0 O < p(x) < A)}, then the mapping n : [-0o,0o] xffl^ 
T> D (f2) defined by IL(p,£) = £ - pKO is a C 2 diffeomorphism. Moreover Dn(0^)(l,e) = 
e — i>£ for any e belonging to the tangent space Tj(90) to dil at £. If x G ^(fi), we 
define the reflection of x through dVl by ip(x) = £ x + p( x ) 1/ i a - Clearly ip is an involutive 
diffeomorphism from £1 fl T 1 ^,, (f2) to O c n TJg^fi), and Dif)(x) = I for any a; G 90. If a 
function w is defined in Q, n T^ (J7), we define w in Tjg^O) by 

( ifarefinT^fi) 
■u(x) = < (2.4) 

( -co^i) if x g fi c nT ft (fl). 



Lemma 2.4 Assume that G 90. Let v G C 1,a! (fi n Tp (Q) \ {0}) be a solution of fTT")) in 
fl (~l T]g (f2) vanishing on dfl \ {0}. TTien S G C ,a (TJg(f2) \ {0}) is solution of a quasilinear 
equation 

Y,^ A ^Dv) = (2-5) 



in Tp(n) \ {0} where the Aj are C 1 functions defined in Tp(Cl) where they verify 

( (i) 1^,0) =0 



(2.6) 



(Hi) 



(x,rf) 



<T\ V 



p-2 



for all x G Tp(n) \ {0} for some [3 G (0,/? ], r\ G R , £ G M and some T > 
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Proof. The assumptions <|2.6 I) implies that weak solutions of (|2.5 I) are C 1,Q , for some a > 
[T^ and satisfy the standard a priori estimates. As it is defined the function v is clearly C 1 in 
7> o (ft)\{0}. Writing Dv(x) = -D(voiP(x)) = -Dip(x)(Dv(ip(x))) and x = ip(x) = ip~ l {x) 

\Dv\ p ~ 2 Dv.DCdx 

\DiJj(Dd)\ p ~ 2 Di/j(Dv).Dip(D() \Dip\ dx. 



in n7>(n) 

But 

^dtpi dv \ I s-^dipj d( 



D^(Dv).Di>(D() = J2 E 



dxk dxi } \ ^-r 1 dxk dxj 
/ \ j 



3 

> dipi dipj dv \ d( 
1 dxk dxk dxi j dxj 

i,k I J 



We set b(x) = \Dip\, 



and 

A(x,r?) - (A^ry), .., A* for?)) - |£ty| I^WI"" 3 (2.8) 

For any £ G <9S1, the mapping -DV'SfKC) is the symmetry with respect to the hyperplane 
tangent to <9S1 at £, so D-i/KOI = 1- Inasmuch is continuous, a lengthy but 
standard computation leads to the existence of some G (0, /3q] such that (|2.6 I) holds in 
7>(si)nsi . if we define A to be \rj\ p 2 77 on T / g(f2)nJl and A on Tp(fl)C\£l C , then inequalities 
HHP are satisfied in T (Q). □ 

These three results allows us to prove our main result 

Theorem 2.5 Let S7 be an open subset ofM. N with a compact C 2 boundary, p(x) = dist (x, <9S1) 
and a G <9S1. Then there exists one and only one positive N -harmonic function u in SI, van- 
ishing on dfl \ {a} verifying 

lim \x — a\ u(x) = — (a, n a ) (2-9) 

x — > a 
]§G§T -> " 

uniformly on S' A '~ 1 n SI, and 

= o(ln|x|)) as x — > oo, (2-10) 

if fl is not bounded. 

Proof. Uniqueness follows from l|2.9 [) by the same technique as in the previous propositions. 

— i 2 

Step 1 (Existence) . If SI is not bounded, we perform an inversion X m with center some 

i — i 2 i — i 2 

to G f2. Because of (|2.10 L the new function u oI„ a| is TV-harmonic in SI' = I™ a (SI) 

and satisfies (|2.9 () . Thus we are reduced to the case were SI is bounded. Since SI is C 2 , 

it satisfies the interior and exterior sphere condition at a. By dilating SI, we can assume 

that the exterior and interior tangent spheres at a have radius 1. We denote them by 
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Bi{uj e ) and Bi(u! z ), their respective centers being uf = a — n a and u e = a + n a . We set 
V l (x) = U l (x — lu 1 ) and V e (x) — U e (x — uo e ) where U % and U e are the two singular N- 
harmonic functions described in Proposition 12 . 21 and Proposition ^. 31 respectively in Bi(lu 1 ) 
and B\ (w e ), with singularity at point a. For e > 0, we put f2 e = 0\ B e (a), S e = fl dB e (a) 
and 9*57 e = <9f2 (~1 -B^r (a). Let u e be the solution of 

diw(|£'u e | A ' _2 Du e ) = in f2 e 

u e = on d*il e (2.11) 
u e = V e on E £ . 

This solution is obtained classicaly by minimisation of a convex functional over a class of 
functions with prescribed boudary value on <9f2 e . For any x G Bi(cu l ), there holds 

dist [x,dB x (u e )) = \x - uo e \ - 1 > dist (x,dfl) > dist (x, dB^uj 1 )) = 1 - . 

thus 

V^ ; (x) < T/ e (a;) VseB^w'), 
by using l|2.1 |l . (|2.3 I) and the maximum principle. Therefore 

V l {x) < u e (x) < V e {x) \fx G n fi £ 

and 

< V e (x) Vx G fi e . 
Finally, for < e' < e, u e ,\ Se < V% e = u £ | Se . Thus 

u t i{x) < u e (x) Vi £ O f . 

The sequence {u e } is increasing with e. By classical a priori estimates concerning quasilinear 
equations, it converges to some positive iV-harmonic function u in Q which vanishes on 
dtt \ {a} and verifies 

V l (x)<u(x) VieBi(w'), 

and 

u(x) < u e (x) Vx g n. 

This implies 

l—t^L<u(x) VigBjK), (2.12) 



2 b-a 



\x - u! P \ 2 - 1 



"0*0 < J 2 Vz G fi, (2.13) 

2 |ir — a| 

By scaling we can prove the following estimate 

u (x) < C P ^ X \ Vx G fi. (2.14) 
|x — a| 

for some C > 0: for simplicity we can assume that a is the origin of coordinates and, for 
r > set u r (y) = u(ry). Clearly u r is ./V-harmonic in Q/r and 

max{\Du r (y)\ : y G fi/rn(B 3/2 \B 2/3 )} < C*max{|u r (z)| : z G fi/rn(B 2 \B 1/2 )}, 

where C, which depends on the curvature of dQ/r, remains bounded as long as r < 1. Since 
Du r (y) = rDu{ry), we obtain by taking ry = x, \y\ = 1 and using (|2.13 \ with general a, 
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|Du(a:)| < C \x — a|~ 2 . By the mean value theorem, since u vanishes on dfl \ {a}, i|2.14 I) 
holds. 

Step 2. In order to give a simple proof of the estimate l|2.9 \ . we fix the origin of coordinates 
at a — and the normal outward unit vector at a to be — e^v- If u is the extension of u 
by reflection through dil, it statisfies ll in Tp(il)\{0} (see lemma 12.4 (I . For r > 0, set 
u r (x) = ru(rx). Then u r is solution of 




(2.15) 



in Tp/ r (Q/r)\{0}. By the construction of Aj(x,r]), we can note that 



lim A> (rx,rj) = \ri\ p ~ 2 Vi, e R N . 

r— >0 

Furthermore, for any x € T ( g(ri)\{0}, ^(x) = p(ip(x)) and c|x| < |V'( a; )l < c_1 M f° r some 
c > 0, the estimate (|2.14 II holds if u is replaced by u r , Q by Tp/ r (Q,/r) and p(x) by /O r (x) := 
dist(x, f2/r) i.e. 

|ti r (x) < C|x|- 2 p r (x) Vx e T 0/r {9,/r). 

For < a < 6 fixed and for some < r$ < 1 the spherical shall T a ,b = {x G : a < \x\ < b} 
is included into Tp / r (fi/r) for all < r < ro . By the classical regularity theory for quasilincar 
equations J2| an d lemma 12.4 I there holds 

H^ r |lc=(r 2/3 , 3/2 ) < C r ll« r llL»(r 1/2 , a) . 

where C r remains bounded because r < 1. By Ascoli's theorem, (|2.12 It and 12.14 \ , u r (x) 
converges to xn\x\~ 2 in the C 1 (r 2 /3.3/2)-topology. This implies in particular that r 2 Du(rx) 
converges uniformly in T 2 /z,s/2 to — 2xjv|x|~ 4 x + x|~ 2 e_/v. Using the expression of Du in 
spherical coordinates we obtain 

r 

r 2 u r i — ru^e H V^u — > — 2cr/vi + e^r uniformly on S N ^ X as r — » 0, 

sin 

where cos0 = xat|x| _1 , i = a;/|a;|, e is derived from x/\x\ by a rotation with angle it/2 in 
the plane 0,x,N (N being the North pole), and V CT / is the covariant gradient on S N ~ 2 . 
Inasmuch i, e and V CT ' are orthogonal, the components of e^r are cos </>, sin </> and 0, thus 

ru^(r, cr', <f>) — » — sin0 as r ^ 0. 

Since 

u{r,a',(t>)= I u c i > (r,a / ,9)d9, 

J-k/2 

the previous convergence estimate establishes <|2.9 II . □ 



Definition 2.6 W^e s/iaH denote by u\ ta the unique positive N -harmonic function satisfying 
\2.9 )) . and call it the fundamental solution with a point singularity at a. 
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3 The classification theorem 



In this section we characterize all the positive TV-harmonic functions vanishing on the bound- 
ary of a domain except one point. The next statement is an immediate consequence of 
Theorem EH and Th. 2.11]. 

Theorem 3.1 . Let Q be a bounded domain with a C 2 boundary and a £ dfl. If u is a 
positive N -harmonic function in f2 vanishing on d£l \ {a}, there exists M > such that 

u(x) < Mui ta (x) \fx e fi (3.1) 

In the next theorem, which extends [21 Th. 2.13], we characterize all the signed N- 
harmonic functions with a moderate growth near the singular point. 

Theorem 3.2 . Let Q, be a bounded domain with a C 2 boundary and a € dQ. Assume that 
Ui ta has only a finite number of critical points in fi. If u is a N -harmonic function in 
vanishing on dtt \ {a} verifying \u(x)\ < Mu\ ta {x) for some M > and any x € fi, there 
exists k e [— M, M] such that u — ku\^ a . 

Proof. We define k as the minimum of the t such that u < lu\ a in f2. Without any loss of 
generality we can assume k > 0. Then either the tangency of the graphs of the functions u 
and fcui. Q is achieved in f2 \ {a}, or it is achieved asymptotically at the singular point a. In 
the first case we considered two sub-cases: 

(i) The coincidence set G of u and ku\, a has a connected component uj isolated in CI. In this 
case there exists a smooth domain U such that ZU C U and 5 > such that ku\, a — u > 5 on 
dlA. The maximum principle implies that kui >a — u > 5 in U, a contradiction. 

(ii) In the second sub-case any connected component u> of the coincidence set touches 9ri\{a}, 
or the two graphs admits a tangency point on dtl \ {a}. If m G uj D dQ \ {a} or is such 
a tangency point, the regularity theory implies du(m)/dii m = ku\^ a (m)/dn m . By Hopf 
boundary lemma, ui^ a {m) / dn m < 0. By the mean value theorem, the function w — ku\ ta — u 
satisfies an equation 

Lw = (3.2) 

which is elliptic and non degenerate near m (see 0, ^), it follows that w vanishes in a 
neighborhood of m and the two graphs cannot be tangent only on <9f2 \ {a}. Assuming that 
uj fl, let xq S il \ uj such that dist (xq, uj) = rg < p(xo) = dist (xo, dfl, and let yo € to be 
such that \xq — yo\ — ^o- Since Ui >a has at most a finite number of critical points, we can 
choose xq such that yo is not one of these critical points. By assumption w = ku\^ a — u 
is positive in 73 ro (xo) and vanishes at a boundary point yo- Since the equations are not 
degenerate at yo there holds 

kdui, a (y )/dv - du(y )/du < 

where v = (yo — xo)/ro, which contradicts the fact that the two graphs are tangent at yo- 

Next we are reduced to the case where the graphs of u and fcui )0 are separated in Q and 
asymptotically tangent at the singular point a. There exists a sequence {£ n } C such 
that limn^oo M(£ n )A*i, (£n) = k. We set \x n - a\ = r n , u n (y) = r n u(a + r n y) and v n (y) = 
r„wi,„(tt + r n y). Both u n and v n are TV-harmonic in f2 n = (f2 — a)/r n . The functions u n 
and v„ are locally uniformly bounded in f2„ \ {0}. It follows, by using classical regularity 
results, that, there exists sub-sequences, such that {w n J and {v nk } converge respectively 
to U and V in the Cj^-topology of fl Uk \ {0}. The functions U and V are A^-harmonic 
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in H « = {x = {x\,X2, •••,£Ejv) : > 0} and vanish on dH \ {0}. Since it can be 
assumed that (£ nfc — a)/r 1lk — > £, there holds [/ < fcV in if, £/(£) = A;V(£), if £ G if, and 
dU(£)/dx N = kdV{C)/dx N > 0, if £ G <9ff (notice that |£| = 1). If £ G 9ff, Hopf lemma 
applies to at £ and, using the same linearization with the linear operator L as in the 
previous proof, it yields to U — kV. If £ G if, we use the fact that \Dux >a (x)\ > (3 > 
for | a; — a| < a for some /3, a > 0. Thus |f>u n (£)| > (3. The non-degeneracy of V and the 
strong maximum principle lead again to U = kV. Whatever is the position of £, the equality 
between U and kV and the convergence in C\ oc leads to the fact that for any e > there 
exists n e G N such that n > n e implies 

(k — e)ui j0 (x) < u(x) < (k + e)ux ta (x) Vx G fl D dB Tn (a). 

By the comparison principle between iV-harmonic functions this inequality holds true in 
fi \ 9-B r „ (a) ■ Since r n — > and e is arbitrary, this ends the proof. □ 

Remark. The assumption that u± ta has only isolated critical points in fl is clearly satisfied 
in the case of a ball, a half-space or the complementary of a ball where no critical point 
exists. It is likely that this assumption always holds but we cannot prove it. However the 
Hopf maximum principle for p-harmonic functions (see implies that U\, a cannot have 
local extremum in Q. 



4 Separable solutions of the ]>harmonic spectral prob- 
lem 

In this section we present a technique for constructing signed TV-harmonic functions, regular 
or singular, as a product of functions depending only on one real variable. Some of the 
results were sketched in |T5]- The starting point is the result of Krol [S] dealing with 
the existence of 2-dimensional separable p-harmonic functions (the construction of singular 
separable p-harmonic functions was performed in jl]). 

Theorem 4.1 (Krol) Let p > 1. For any positive integer k there exists a unique (3k > 
and iOk '■ K i— > M, with least antiperiod w/k, of class C°° such that 

u k {x) = \xf k u k {x/\x\) (4.1) 

is p-harmonic in M 2 ; (3k is the unique root > 1 of 

s , pk 2 + (p - 2)(2k - 1) , , N 

(2k-l)X 2 -- — - -X + k 2 = 0. (4.2) 

p- 1 

((3k, Wk) is unique up to translation and homothety over ujk- 

This result is obtained by solving the homogeneous differential equation satisfied by ujk = u>: 

' 2 ^ 2)/2 *e) = & (1 + 09 - l)(p - 1)) (/? 2 - 2 + c e 2 ) (p - 2)/2 w . (4.3) 

In the particular case k = 1, then (3\ = \ and cl>i(0) = sin 6*. For the other values of k the 
(3k are algebraic numbers and the u)k are not trigonometric functions, except if p = 2. More 
generally, if one looks for p-harmonic functions in M. N \ {0} under the form u{x) — u(r, o~) — 
r^v(a), r = \x\ > 0, a = xj \x\ G S N ~ X , one obtains that v verifies 



-div r 



( (/3V + |V CT w| 2 J V CT t> j = Ajv,/3 (/?V + I V ctW | 2 J v (4.4) 
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on S N ~ X , where Ajv,/3 = ft (N — 1 + (/3 — l)(p — 1)) and rfw ff and V CT are respectively the 
divergence and the gradient operators on S N_1 (endowed with the Riemaniann structure 
induced by the imbedding of the sphere into R w ). This equation, called the spherical 
p-harmonic spectral problem, is the natural generalization of the spectral problem of the 
Laplace-Beltrami operator on S*^ -1 . Since it does not correspond to a variational form 
(except if p = 2), it is difficult to obtain solutions. In the range oil < p < N— 1, Krol proved 
in the existence of solutions of 1)4.4 [I . not on the whole sphere, but on a spherical cap 
(which reduced (|4.4 I) to an non- autonomous nonlinear second order differential equation). 
His methods combined ODE estimates and shooting arguments. Later on, Tolksdorf [TT] 
introduced an entirely new method for proving the existence of solutions on any C 2 spherical 
domain S, with Dirichlet boundary conditions. Only the case (3 > was treated in 
and, by a small adaptation of Tolksdorf approach, the case f3 > was considered in JSj ■ We 
develop below a method which allows to express solutions as product of explicit one variable 
functions. 



4.1 The 3-D case 

Let (r, 9, 4>) € (0, oo) x [0, 2ir] x [0, tt] be the spherical coordinates in 



x\ = r cos a sin < 
x 2 = r sin 9 sin (, 
x 3 — r cos <p 



Then l|4.4 \ turns into 



2 \ (P"2)/2' 



sin (j) (3 v + vi + 



sin 



We look for a function v under the form 

u(M) 



then 



/?V +vl + 



sin 



[3 2 v 2 + vl + 



/3(2+(/3-l)(p-l)) [P 2 v 2 +vl + - 2 

sin c 



2 \ (p-2)/2 - 



t'0 



sin 

2 \ (P-2)/2 



(sin^w^) 
(sin0) 2 ^ 2 (/3 2 ^ 2 + c, 2 ), 



(4.5) 
(4.6) 



sm 



and 



S11K 



2 \ (p-2)/2 



v sm 



o2 2 , 2 , 
P W + «<A + 



sm 



XI = (sm0) (/3 - 1)(p - 1)+1 (^ 2 + W 2 )^ 2 )/ 2 ^, 



„2 \ (P"2)/2 



2 \ (p-2)/2' 



sin0 /3V +W 2 + —4- W 



: ^(sin^)^- 1 )^" 1 )- 1 [((/? - l)(p - 1) + 1) - sin 2 <j> ((/? - l)(p - 1) + 2)] (/3 2 co 2 + w 2 )(p- 2 )/ 2 w. 
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It follows that lo satisfies the same equation 1)4.3 \ . The next result follows immediately 
from Theorem 14. II 

Theorem 4.2 Assume N = 3 and p > 1. Then for any positive integer k there exists a 
p-harmonic function u in R 3 under the form 

u(x) = u(r,6,<p) = r h (sin <j)f k u k {6) (4.7) 

where (3k and u>k ar & as i n Theorem \4-l\ 

In the case p — 3 we can use the conformal invariance of the 3-harmonic equation in R 3 to 
derive 

Theorem 4.3 Assume p = N = 3. Then for any positive integer k there exists a p- 
harmonic function u m R 3 \ {0} under the form 

u(x) =u(r,6,(f)) = r" /3fc (sin0) /3fc cj fe (6») (4.8) 

where (3k and u>k are as in Theorem \4-l\ with p = 3. 

As a consequence of Theorem 14 . 31 we obtain signed 3-harmonic functions under the form 
()4.7 I) in the half space R 3 ^ = {x : X2 > 0}, vanishing on \ {0}, with a singularity at 
x = 0. They correspond to even integers k. The extension to general smooth domains fHs a 
deep chalenge. In the particular case k — 1, we have seen that (3i = 1 and ioi{8) = sin 6 = X2, 
that we already know. 

4.2 The general case 

We assume that N > 3 and write the spherical coordinates in under the form 

x = {{r,a) E (0,oo) x S^ 1 = (r, singer', cos0) : a' E S N ~ 2 , <p E [0,tt]} . (4.9) 
The main result concerning separable p-harmonic functions is the following. 

Theorem 4.4 Let N > 3 and p > 1. For any positive integer k there exists p-harmonic 
functions in M. N under the form 

u(x) = u(r,a',<j>) = (rsmcf)) 0k w(a'). (4.10) 

where (3k is the unique root > 1 of \4-2 | ) and w is solution of \4-15 | ) with (3 — (3k- Further- 
more, if p = N there exists a singular N -harmonic function under the form 

u(x) =u(r,o-',ct>) = r _,3fc (sin 0)' 3fe w(a'). (4.11) 

Proof. We first recall (see jEj for details) that the SO(N) invariant unit measure on S N ^ X 
is da = on sin ~ 2 (f>da' for some etjv > 0, and 

V a v = -v^e + -?—y a iv, 
sin<p 

where e is derived from x/ \x\ by the rotation of center angle 7r/2 in the plane going thru 
0, xj \x\ and the north pole. The weak formulation of 1)4.4 \ expresses as 

Jo Js N - 2 \ sm 4> J \ sin / 

= ^n,J I (p 2 v 2 +v 2 + -^-\V a ,v\ 2 ) {P ] vCsm N - 2 cj>da'dct> 
Jo Js N -A sm 4> ) 

(4.12) 
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or, equivalcntly 



1 



sin 



N-2 



sin 



N-2 i / fl2„2 , „,2 



1 



sin 



(p-2)/2 



-div a > 



sin 



/3V + t>l + - r - 5 -|V < ,/t>| 



sin 



(p-2)/2 



(p-2)/2 



(4.13) 



Aiv,/3 L8V + «5 + T-2- |V CT ^|' 



sin 



where ditv is the divergence operator acting on vector fields on S N 2 . We look again for 
p-harmonic functions under the form 



u(r, <t) — u(r, a', (j>) — r^v(a', (j>) — r 13 sin' 3 (f)w(a'). 



(4.14) 



Then 



thus 



P 2 v 2 +vl + —^-\V al v 
sin 



, (P"2)/2 

= (sin0)^- 1 )(f- 2 )(/3 2 U ; 2 + |V o 



(p-2)/2 



sin 



iV-2 



sin"" 2 /3V + i; 2 



1 



sin 



|V^«| : 



(p-2)/2 



^(sin^^- 1 )^" 1 )- 1 ((iV — 2 + (/3 — l)(p - 1)) - (JV - 1 + (/? - l)(p - 1)) sin 2 , 

, (p-2)/2 



(p 2 w 2 + |V 
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and 



-div„i 



sin 



r^+^ + -T-| V -'^l' 
v sin <p 



(P"2)/2 



= (sin^)^- 1 )^- 1 )- 1 ^^ 



>-2)/2 



Finally to satisfies 



— dlVa' 



(/3 2 w 2 + |V t 



2 ,(p-2)/2 



(/3W + |V„H 2 ) 
Xn-i,0 [f3 2 w 2 + \V a >w\ 2 ) w (4.15) 



on 5 , which is the desired induction. 



□ 



In order to be more precise, we can completely represent the preceding solutions by 
introducing the generalized Euler angles in M. N = {x = (x\, ...,xn)} 



x\ = r sin#/v-i sin0jv_ 2 ... sin# 2 sin#i 
X2 = rsin0jv-i sin0jv_ 2 ... sin# 2 cos^i 



xjv_i = r sin f jv-i cos Vn-2 
x N = rcos0N-i 



(4.16) 
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where 9% 6 [0,2tt] and Ok G [0,7r], for k = 2, ...,7V — 1. Notice that 9m~i is the variable <j) 
in the representation l|4.9 I) . The above theorem combined with the induction process yields 
to the following. 

Theorem 4.5 Let N > 3 and p > 1. For any positive integer k there exists p-harmonic 
functions in M. N under the form 

u(x) = (rsm0 N _ 1 sm0 N _2...sm0 2 )' 3k uJk(9i) (4-17) 

where ((3k,0Jk) are obtained in Theorem \4-l\ Furthermore, if p — N there exists a singular 
N -harmonic function under the form 

u(x) =r-P k {sm6 N - 1 sm6 N - 2 -si-n&2)' 3k Uk(0i). (4.18) 
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